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ABSTRACT 

The molecular cloud, DR21 Main, is an example of a large-scale gravitational 
collapse about an axis near the plane of the sky where the collapse is free of major 
disturbances due to rotation or other effects. Using flux maps, polarimetric maps, 
and measurements of the field inclination by comparing the line widths of ion 
and neutral species, we estimate the temperature, mass, magnetic field, and 
the turbulent kinetic, mean magnetic, and gravitational potential energies, and 
present a 3D model of the cloud and magnetic field. 

Subject headings: Molecular Clouds, Submillimeter, Polarization, Star Formation 



1. Introduction 

Partially ionized material in a molecular cloud will respond to gravitational forces by 
moving, primarily, along magnetic field lines. If one assumes that the field is not significantly 
tangled by turbulence, one can expect a concentration of material in an obla te core. This 



mode l ha s been verified by t he analysis of polarimetric and photometric data (IDotson et al. 



20091 ) by iTassis et al.l (120081 ) showing that the field is directed along the short axes of the 
cores. The neutral component of the material will move slowly across field lines toward 
the center of gravity thus increasing the gravitational field. Where a pinch in the field is 
indicated by a polarization map, one can infer that the gravitational field in the core is 
competing with the support provided by the field. Within a critical radius, one expects the 
cloud to collapse. 



^Dissertation submitted to the Department of Astronomy and Astrophysics, University of Chicago, in 
partial fulfillment of the requirements for the Ph.D. degree. 
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Examples of p inched fields have been found in Orion fjSchleuninpl998l ) and in NGC1333 
( IGirart et al.ll2006l ). In the case of DR21 Main (DR21M), the polarization map extends to 
radii such that the field straightens out as the gravitational field diminishes leaving an 
hourglass field configuration. In this cloud one has an opportunity to estimate the mass 
within the collapsing region. We use 350 /im polarimetry to measure the field as projected 
on the sky ( §3.11) : line width measurements to estimate the inclination of the field to the line 
of sight ( §4.11) : photometric maps to estimate the mass distribution ( §4.3p : and the angular 
dispersion of the field vectors to estimate the field strength ( §4.41) . 



2. DR21 



The giant star formi ng complex, DR21 (Fi gure [T]) is located in the Cygnus constella- 
tion ~ 3kpc from Earth (ICampbell et al.lll982l). The southernmost component, DR21M, 



has a mass of ~ 20, 000 Mq (IRichardson et al.l Il989l). It contains one of the most en 



ergetic star f o rmati on outflows detected (I Garden et al.l Il99ll : iGarden fc CarlstromI Il992l ) . 
Garden et al.l (Il99ll ) measured the mass of the outflow to be 3000 Mq. Observers using 
Spitzer fou nd 5 near-infrared (NIR) sources (see Figure [2]) inside DR21M, presumed to be 
protostars (IKumar et al.l 120071 : iDavis et al.l 120071) indicating that regioi is within the cloud 
are un dergoing gravitational collapse. iGlenn et al.l (119991 ) at 1.3 mm and lMinchin fc Murray 
( Il994j ) at 800 /im have previously measured the polarization of DR21 finding uniform po- 
sition angles ~ 25° East of North (i.e. B field ^ 65 "W est of North). The galactic plane 



is oriented at 42 °East of North. iRoberts et al.l (119971 ) measured a line-of-sight magnetic 
field of ~ 400yuG using HI Zeeman measure ments around t he du st emission peak and star 
formation regions. The regions observed by iRoberts et al. (|l997h with 3.5(7 detections are 
outlined in green in Figure [31 



3. Observations 



3.1. Polarization 



Obse rvations, made at the Caltech Submillimete r Observatory (CSO) with the polarime- 
ter Hertz (ISchleuning et al.l 1 199 71 : iDowell et al.lll998l). provided 3 50 /im polarization maps of 



the region. Hertz was later re placed by SHARP (ILi et al.ll2008l ). a module mounted ahead 
of the photometer SHARC-II (jDowell et al.ll2003l ) to permit polarization measurements. 



SHARP has a wire cross-grid to split the incoming signal into orthogonal components 
of polarization which are then detected on two 144 pixel subsets of the SHARC-II detector 
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Fig. 1. — Hertz image (jDotson et al.l 120091 ) of the star forming region DR21. Red vec- 
tors (E-vectors) are for measurements better than 3a polarization and blue vectors are for 
measurements of 2a to 3a. Open circles are locations with a 2a upper limit of 1% 
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Fig. 2. — SHARCI I 350 /xm flux ma p with x's marking locations of Spitzer NIR sources 
jKumar et aDbnO?! : [Pavis et aliboOTh 
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array. Standard infrared techniques of chopping the secondary mirror to remove sky emission 
and positiona l offsets to remove grad ients in sky emission were used for all the polarimetry 
observations (iHildebrand et al.ll2000l ). The chopping direction was azimuthal with a throw 
of 300". To avoid chopping into the adjacent cloud, DR210H, all observations were made 
when the chop throw was > 30°from north-south on the sky A half-wave plate was rotated 
to four different angles separated by 22.5°. Each observation cycle consisted of a set of 
observations at the 4 half-wave plate angles. The signals from the two 144 pixel arrays were 



then combined to give the polarization and photometric signals as described by iKirby et al. 



(120051 ). The polarization and flux signals are given by 



Polarization Signal = H — V (1) 
Flux Signal = H + V (2) 

where H and V are the horizontal and vertical components of the signal. The Stokes param- 
eter, Q, was calculated from the difference of the polarization signals of the 1st and 3rd angle 
in the cycle, while the parameter U was calculated from the difference of the polarization 
signals of the 2nd and 4th angle. The source flux, /, was determined by the average of the 
flux signals from the 4 half- wave plate positions. The Q, f/, and / from each cycle were then 
combined after correcting for changes in attenuation and sky noise. Due to the lack of an 
instrument rotator, this combination had to be done a new way. The corrected Q, U, and 
/ from all of the cycles were combined to form an irregular ly sampled map. This map w as 



then smoothed to a resolution of ~ 10" (for more detail see iHoude fc VaillancourtI (120071 ) 



The observations were made in 7 fields spaced by 45". Each field was observed with 
a 4-point dither on corners of a square 10" wide in Right Ascension and Declination. A 
few additional observations were taken to fill in areas of low signal to nois e in the above 7 



fields. The data were corrected for instrument polarizations of ~ 0.3 — 0.5% (ILi et al.l (120081 ). 
Vaillancourt et al, in preparation) 

The SHARP observations were taken on 2006 July 15-19 during intermediate weather 
(t225GHz ~ 0.07 — 0.08 corresponding to 1.7-2.0 mm zenith water vapor). The SHARP 
polarization data set shows 1270 sky positions with a signal to noise > 3. The beam size 
of SHARP is ~ 4x the pixel size. Limiting the sample so there is no beam overlap leaves 
78 measurements. The median value of the polarization is 2.4% ranging between 0.37% and 
13%. Half of the points have a signal to noise better than 4.7. At this value, the accuracy of 
the position angle of the polarization vectors is ~ 6.5°. The direction of polarization shows a 
pinched or hourglass shape (see Figure [3]). The axis of symmetry of the hourglass is inclined 
~ 15° East of North (determined by model see | 



The results tabulated by 



Dotson et al, 



(12009! ) were obtained with the polarimeter. Hertz 
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Fig. 3. — 350 fim SHARP(red) and Hertz(blue dotted) polarization of DR21M. B vectors are 
all set to the same arbitrary length in the magnetic field direction. The grayscale background 
comes from a SHARC-II fiux map. Gre en contours show the locations larger than a SHARC- 
II pixel with 3.5a Zeeman results from [Roberts et al.l (119971 ) 
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(Figure[3]blue dotted vectors). The Hertz measurements were taken on three observing runs, 
1997 April 18-27, 1997 September 18-26, and 2001 April 10-16. 



3.2. Heterodyne Measurements 

The polarimetry measurements from SHARP and Hertz are supplemented with 345 
GHz heterodyne measurements of ion and neutral molecular line widths from the CSO. The 
heterodyne measurements were taken on 2006 July 12-13 {t225ghz ~ 0.1). The inclination of 
the field at eleven points within the pinched region was determined from measurements of line 
widths of neutral and ionized molecules of comparable mass. The heterodyne measurements 
of the J— »4-3 transiti on in HON and HCO-I- were made using only on/off position differencing 
jHoude et all boOOal lbh. The data were reduced using the CLASS program as part of the 



GILDAS package of programs (http://www.iram.fr/IRAMFR/GILDAS). Each of the lines 
was smoothed to half of the natural resolution using the CLASS routine SMOOTH to improve 
the signal to noise. One of the positions was removed due to low signal to noise. Another 
position was unusable because bandwidth of the instrument was insufficient to image the 
entire line. 

The spectrum at each position was then fit by a multi-Gaussian model (Figure HI) to 
determine the line widths (more properly the standard deviation a^). Multiple Gaussians 
were used to account for the self-absorption seen in most of the positions, and for one- 
sided outfiows, which can cause an artificial reduction in line widths if not accounted for 
properly. The effects from broad outfiows were removed from two locations. The fits 
were done using the XGaussfit program from the FUSE package of software for IDL ( 



|http: / / fuse.pha.jhu.edu/analysis/fuse_idl_tools.html ) . Results of the fits are given in Table 



m 



4. Results 
4.1. Inclination Angle 



Houde et all J2000al lbl) investigated the effect of a magnetic field on the velocities and 



ions in a partially ionized fiow (not necessarily an outfiow or jet but any situation where 
local mean velocity is not zero). For a significantly strong field, the ions were forced into 
gyromagnetic rotation about the magnetic field instead of following the fiow. This rotation 
led to a reduction in line width and suppression of high velocity wings. The line widths 
were derived for both the neutrals and ions. The ratio of the linewidths depends only on the 
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Table 1. Line Width Data and Inclination Angle with the Line of Sight for DR21 
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Fig. 4. — Line profile (black) and best fit multi- Gaussian model (red) of the J^4-3 transition 
of HCN (top) and J^4-3 transition of HCO+ (bottom) at the center of DR21M (20'^37"14.1^ 
42°8'53"). 
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orientati on of the neutral flow(s) and the inchnation angle of the magnetic field with the line 



of sight (iHoude et al.l 120021 . |2004J ). Assuming that the neutral component of the material is 
H2 (mean molecular mass of 2.3), the square of the ratio of the ion to neutral line width is 
given in terms of the angle with the line of sight, a, by 

afi ^ e cosa^ + /(O.lGcosa^ + 0.84sinaV2)(^i/Aii - 1)"^ 

~ (3) 



af,^ e cosa^ + /sinQ;^/2 



where 



cos^Ae 

(4) 



6 

, 2-3cosA0 + cos3Ae 

/= -Q , (5) 

and rrii and fii are the ion mass and the reduced mass respectively. The neutral flows are 
modeled as symmetrical and contained within a cone of width centered on magnetic 
field direction (symmetry axis, see Figure [5]). All values of a must satisfy t he condition 



P/ Prn nT. < slu wherc Pmax IS the maximum polarization seen by Hertz (10 % lHoude et al. 



(12004)). 



Although the cone half width, A^, can, in principle, be determined by finding the curve 
that best fits the above condition on P/P^ax^ the data for DR21 do not provide significant 
constraints on this quantity. The values of a presented in Table [1] are computed assuming 
A^ = 45°. The inclination with the line of sight decreases toward the center of DR21 as 
would be expected for a magnetic field in the center of a gravitationally contracting region 
with an axis near the plane of the sky. Inclination angles for points located at the edge of 
the cloud, where the gravitational pinch has had only a small effect, should approximate the 
inclination of the cloud's mean B field. These points for DR21 give an estimate of ~ 70° 
inclined to the line of sight or ~ 20°to the plane of the sky. 

Choosing the wrong value of A6' leads to a systematic error in the calculation of each 
angle, but the difference between angles is not significantly affected. The overall trend 
of decreasing angle toward the center of the cloud remains. If, instead, we had assumed 
AO = 90°, the values shown in the table would decrease by ~ 15° but the overall trend 
toward decreasing inclinations toward the center would not change. 



4.2. Temperature and Optical Depth 



The spectral energy distribution at 350 and 850 jjm was fitted to give estimates of 
the dust temperature and optical depth at 350 /im. The 850 /im map was acquired from 



Fig. 5. — Illustration of flow of neutral material. All flows (velocity vectors) are contained 
within a co ne of width cent ered on the symmetry axis deflned by the magnetic fleld (see 
Figure 2 of iHoude et al.l ((2002) • Length of vectors is arbitrary as magnitude and direction 
are assumed independent. 
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Fig. 6. — Polarization vectors with known 3- dimensional spatial orientation. The angle of 
the vectors is drawn at the angle of the magnetic field in the plane of the sky. The length 
of the vector is proportional to the angle, a, of the magnetic field with respect to the line of 
sight. The grayscale background is from SHARC-II. 
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the SCUBA online archive @ (Project ID: m02bu47). The SHARC-II map at 350 fim was 
smoothed and repixelated to match the 850 /im map. The intensity at frequency u and 
corresponding wavelength A was then modeled as 

J, = B,(T)(l-e-"35°(=^^°/^)'), (6) 

where the exponent gives the atmospheric transmission at frequency u as determined from 
the transmission at 350 /xm, and B,j{T) is the Planck function at frequency v and temperature 
T. An example fit to the data for the center of DR21 is shown in Figure [3 This formula 
assumes the cloud is optically thin at the wavelengths used. The ratio of the the 350 /im 
flux to the 850 /im is constant throughout the majority of the cloud (~ 15) with a discrete 
jump at the edges of the cloud (~ 25) and a dip in the center (~ 8). The constancy of the 
ratio implies optical thinness except at the center where the ratio drops. 

Temperatures at the peak were ~ 20 K and went up to ~ 25 K in the cloud. The edges 
of the cloud where the ratio increased had higher temperature (30-50 K). The 25 K region of 
the cloud extended along a rough northeast-southwest direction similar to the direction of the 
outflow. The optical depth showed a similar trend with a value at the peak greater than 1. 
The optical depth quickly drops below 1 indicating that all but the brig ht peak areas o f DR21 



main are optically thin at 350 /xm, as expected from the flux ratios. iHarvey et al.l (119861 ) 
found temperatures from 50 /im and 100 /zm flux maps to be ~ 50 K. These wavelengths 
may be coming from a differe nt environment w i th a d ifferent temperature. Supporting a 



multi-environment hypothesis, [Richardson et al.l (119891 ) found a 350 /xm to 800 /xm ratio of 



16 at DR210H with a corresponding temperature of ~ 25 K in regions given to be ~ 37 



K by lHarvey et al.l (Il986l ). Throughout the rest of the paper, we used temperatures of 20-25 



K from the SHARP and SCUBA fluxes. 



4.3. Mass 



The mass of cloud within a specified region o f a cloud can be determined from its 
temperature and far infrared flux (IHildebrandl 1 1 9831 ) 



M 



BuiT) 



M 



(7) 



where F is the flux, D is the distance to the cloud, and 

Cm = [N{H + H2)/T{u)]mHli = 1.2 ■ 10^^ {350 /m)'^mHf^, 



^http://www2.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/jcnit 
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where r(z/) is the optical depth of the cloud at frequency v. Within a radius of ~ 1 pc (size 
of SHARP map in E-W direct ion), DR21 main has a m ass of ~ 25, OOOM©, comparing nicely 
to the ~ 20, OOOM© found bv iRichardson et al.l Jl989h . 



4.4. Magnetic Field Strength 



Chandrasekhar fc Fermil (Il953l ) described a method to estimate the magnetic field in the 
plane of the sky. This method related the dispersion of the direction of starlight polarization 
from a straight line to the strength of the magnetic field, 



-^planeofsky \ „ (9) 
V 3 (J/3 

where p is the density of the gas, crobs{v) is the observed line of sight velocity dispersion, and 
(T/3 is the dispersion in polarization direction in the plane of the sky. The weaker the field, 
the greater the dispersion of the polarization vector. 

The dispersion of the polarization vectors will be overestimated if it is taken with respect 
to the mean field. To determine the dispersion about the local magnetic field one must take 
into account the inclination of the mean field ( §4.11) and the structure of the large scale 
field due to non-turbulent effects such as gravitational collapse, expanding HIT regions, and 
differential rotation (Hildebrand et al. submitted to ApJ). A polarization map usually shows 
a smoothly-varying pattern of vectors, therefore at separations small compared to the cloud 
diameter, the 2-point angular correlation function (a/<A/5^ >) of position angles should 
increase almost linearly as shown in the bottom panel of Figure [81 The linear portion would 
have a zero-intercept if there were no measurement error and no turbulence. Both effects will 
cause the whole 2-point correlation to be displaced upward (except at values < correlation 
length of turbulence which is expected to be much less than the SHARP or Hertz resolution 
of DR21). The y-intercept, b, of the best fit line is then the quadratic sum of the estimate of 
the dispersion, (j/3 and the measurement error, (7^6. A cross correlation of repeat observations 
gives the estimate of the measurement error (value at separation) and an estimate of the 
quadratic sum of measurement error and turbulent dispersion (intercept of best fit line to 
linear region at small non-zero separation) 



4 + (10) 



Figure [HI shows the correlation function for two subsets of the Hertz data set. The subsets 
were constructed by taking every other raw data file starting with file 1 for subset 1 and file 2 
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for subset 2. The line best fit to the hnear region of the correlation function has an intercept 
of 9?3. The value of the correlation function at zero separation is 7?6 giving an estimate of 
5?3 for the turbulent d ispersion. Usiri g the C handrasekhar and Fermi method formula with 
correction factor from lOstriker et al.l (120011 ) gives a plane of the sky magnetic field of 3.1 
mG where aobsiv) is the observed line of sight velocity dispersion of HCN measured to be 
4.2kms~^ and p is the mean density calculated from the mass in §4.31 to be 4.4xl0~^^gcm~^. 

The method of data acquisition of SHARP prevents repeat observations of the same 
location (caused by lack of instrument rotator) preventing that sample from being split as 
done for Hertz. However, the mean of the individual errors on the angles give a good estimate 
for the measurement error (7?2 to 7?6 for Hertz). Using that method of estimating the mean 
error and fitting a line to the linear region of the correlation function gives a magnetic field 
estimate of 2.5 mG using SHARP data (dispersion of 6?6). 



Using the measurement of 400 /iG for the line of sight magnetic field by [Roberts et al. 



(119971 ) and taking a ~ 3 mG field in the plane of the sky gives a total magnetic field of 
~ 3 mG inclined to the plane of the sky at ~ 10°. This estimate of the angle is probably 
low due to the Zeeman measurement being from HI, which exists at lower density than the 
dust we o bserve in polarizati o n so the field strength is also likely lower there. Taking the 
scaling of iHeiles fc Crutcher I (120051 ) and the density of the region observed in HI as > 10^, 
the angle with respec t to the plane of the sky is < 53°. For a density of lO^cm"^, the 
angle would be ~ 20° [Roberts et al.l (119971 ) assumed that the magnetic field was less than 
1 mG because unless the field was very uniform, the field would have bent somewhere such 
that they detected a 1 mG field. The dust polarization measurements show a field that is 
very uniform implying the field could be larger than 1 mG and not detected in the Zeeman 
observations. 



4.5. 3-Dimensional Model 

A three dimensional model of the fiux and magnetic field was constructed and fit to the 
observed fiux and polarization data. The fiux was modeled as a Gaussian ellipsis of the form 

2;^ -|- 

Fccexpl — — ], (11) 

where x is coordinate along the symmetry axis of the magnetic field, y and z are the other 
two Cartesian coordinates, and a and b are widths of the Gaussians to be fit. The magnetic 
field was modeled as 

(d cx'^ y d ox^ z \ 

^""'Ix^^^^d^ ^y2j^Z^'l^^^£^d^^f^^ j ' ^^^^ 
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where 

/oo 
Fdx. (13) 

This form was chosen to give a magnetic field with the shape of an hourglass that straightened 
at large distance from the center. The magnetic field strength in the x direction in the 
midplane was set proportional to the column density. Assuming flux freezing, the magnetic 
field strength in the x direction for locations outside the midplane was set to a value in the 
midplane that is found by tracing the field line back to the midplane. The field in the y-z 
plane was assumed radial and the strength was found from the strength in the x direction 
and the overall magnetic field direction. The fiux and magnetic field strengths were then 
rotated and summed along the line of sight. The values of a, b, c, and d were then set to 
give the best fit to the fiux and the magnetic field direction determined by the polarization 
measurements. The best fit values were a = 5,b = 15, c = 10, d = 90 (units for a, b, c, d 
are SHARP pixels = 2'.'375, 20°for the angle of the magnetic field with the plane of the sky, 
and 15°East of North for the angle of the projected minor axis on the plane of the sky. The 
fit for the magnetic field direction was quite good (chi-squared of ~ 2). The fiux fit was 
poor as DR21-Main is not an exact ellipse. Figure [9] shows the model vectors along with the 
SHARP vectors. The angle with the plane of the sky is somewhat higher than found in §4.41 
but agrees well with the measured inclination angles farthest from the center (see §4.ip . 



4.6. Gravitational and Magnetic Potential Energies 



Using the mass and magnetic field estimates discussed in §4.3[ §4.4[ and the 3D model 
presented in §4.5^ one can estimate the dependence on radius of the gravitational and mag- 
netic potential energies. The gravitational potential energy of a fiattened centrally condensed 
spheroid at a radius R centered on the peak of the cloud is 



E, 



grav 



3aGM^ 
5R ' 



and the magnetic energy of the uniform field is given by 



E.. 



mag 



6 



(14) 



(15) 



where the values a = 1.2 and b = 0.3 and are calculate d by setting the e nergies equal 
when the mass to magnetic fiux ratio is the critical value (IMcKee et al.lll993l ). B is scaled 
according the model in §4.51 and normalized such that the energy inside 1 pc is the same as 
for the mean field calculated in §4.4[ By com paring large scale p olarization maps to with 
magnetohydrodynamic turbulence simulations, iNovak et al.l ( 120071 ) found that the ratios of 
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the uniform to fluctuating magnetic field in the clouds NGC 6334 and G333.6-0.2 were in the 
range 0.6-2.0. Figure [TO] shows the gravitational and magnetic potential energies (assuming a 
ratio between uniform to fluctuating field of 1) versus radius. Extrapolating the gravitational 
potential energy curve to find the radius at which the two energies are equal gives an estimate 
of ~ 1.2 parsec. The dotted line in Figure [TOl shows turbulent energy calculated from the 
average HON line widths of 4.2 km/s and assuming the velocities are isotropic. The sum of 
these two energies is the dotted-dashed line and intersects the gravitational energy curve at 
0.7 parsecs. One would then expect the critical radius to be between 0.7 and 1.2 parsecs. The 
critical radius should not be regarded as the radius which is suppor ted against collapse since 



the virial theore m can not be applied in this way to magnetic clouds (lMouschoviaslll987l . Il991 



Dib et al.l 120071 ). Figure [T2] shows a plot of the mass to ma gnetic flux ratio versus r adius . 



The ratio is in units of the critical value of the ratio defined by lMouschovias fc Spitzerl (119761 ) 



(16) 



where c$ is taken to be 0.12 (ITomisaka et al.lll988l ). The ratio dips below the critical ratio 
at ~ 1.2 parsec. The mass inside this region is ~ 25, OOOM© ( §4.31) which is approximately 
the entire mass of the cloud. Figure [11] shows the polarization vectors with circles at the 
critical radii with and without turbulent energy. The magnetic field configuration changes 
at around these radii from the pinched structure to the east-west mean direction (see Figure 
[3]). This geometry implies that the cloud is relatively stable against collapse at these radii, 
implying that the magnetic field is providing a significant amount of support against collapse. 
However the presence of protostars are radii greater than 1.2 parsec in Figure [TT] suggests 
that the true picture may be more complex. 



5. Summary 

We present a three-dimensional model of the molecular cloud DR21 (Main). Polarimetry 
at 350 /im has provided a map of the magnetic field as projected on the plane of the sky; line 
observations of ion and neutral molecules have provided measurements of the inclination of 
the field to the line of sight; photometry has provided column densities and 350/im/850/im 
color temperatures. 

A 3-dimensional model of an hourglass configuration magnetic field inclined to the line 
of sight is fit to the observations. We find a pinched field with an axis of 10°-20° from the 
plane of the sky. The mean field strength is ~ 3 mG. The gravitational potential energy is 
equivalent to energies from support mechanisms (magnetic field and turbulence) at a radius 
between ~ O.Spc — l.lpc, a region encompassing ~ 20, OOOM©. 
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Fig. 7. — SED of center region of DR21 with fit from Equation [61 x's denote data points 
and line is fit from Equation [51 
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Fig. 8. — A 2-pt angular correlation function of the polarization angle of DR21 from 2 
separate subsets of Hertz data. The top figure shows the whole cloud and the bottom figure 
show the same function at separations <~ | of the major axis. The black circle shows value 
of the function at separation. 
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Fig. 9. — 350 /iin SHARP polarization map (red vectors) with overplot of projected plane 
of sky vectors from 3D fit. 
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Fig. 10. — The gravitational potential (dashed line and asterisks) is computed from the 
calculation of mass ( §4.31) . The magnetic field energy (solid line) is based on determination of 
magnetic field strength using 2-dimensional CF method and published Ze eman measurement s 
( §4.4p assuming equipartion between the uniform and fluctuating field (INovak et al.l 120071 ). 
Recent results (Hildbrand et al. in preparation) show that equipartition is an overestimate 
and that the ratio of the turbulent to mean field is approximately 10%. The turbulent energy 
(dotted line and crosses) is calculated from the mass and observed HCN line widths. The 
dot-dashed line and squares is the sum of magnetic and turbulent energies. 
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Fig. 11. — DR21M with SHARP magnetic field vectors. Solid circle delineates radius where 
magnetic energy becomes larger than gravitational potential energy. Dotted circle delineates 
radius where the sum of magnetic and turbulent energy becomes larger than gravitational 
potential energy assuming equipartion. The grayscale background is from SHARC-II and 
black 'x"s mark the location of Spitzer point objects. 
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Fig. 12. — Mass to flux ratio in units of the critical ratio between supercritical and subcritical. 
Asterisks are calculations of ratio and dashed line is interpolation/extrapolation 



